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We find that in an opaque quark-gluon plasma, significant particle production at high transverse 
momenta can occur via acceleration of (initially) lower-momentum partons. This new mechanism, 
which is the opposite of parton energy loss (jet quenching), has important implications for observ¬ 
ables, such as elliptic flow, in heavy-ion collisions at RHIC energies and above. 
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Introduction. - There are several indications that 
an opaque, strongly-interacting quark-gluon plasma has 
been created Hii in relativistic heavy-ion collisions 
at the Relativistic Heavy Ion Collider (RHIC). One key 
observable is the large azimuthal momentum anisotropy 
of particle production Sdf, 6j] in noncentral collisions, 
characterized by the elliptic flow coefficient V 2 = (cos 20), 
which requires more than an order of magnitude higher 
opacities than that of a weakly-coupled perturbative par- 
ton plasma Q ■ Though such a dense system is still dissi¬ 
pative 0, i.e., it is not an ideal fluid, transport mean free 
paths A ~ 0.1 fm are as short as their quantum mechan¬ 
ical lower bound |§] and, therefore, the system seems to 
be the most perfect fluid possible. 

The large parton opacities have several consequences. 
For example, charm quarks and antiquarks are predicted 
to acquire a large elliptic flow |n| , and space-momentum 
correlations of energetic particles are expected to show 
a strong surface emission pattern & In this paper 
we show that particle production at high-pr is also af¬ 
fected because there is a considerable probability that a 


low-momentum parton gets accelerated to high transverse 
momenta via multiple scatterings. 

We emphasize that the effect considered here is differ¬ 
ent from the parton energy gain due to detailed balance 
in a static thermal system Il2il . The new source of high- 
Pt particles here is the collective dynamics of the dense 
parton medium. These particles would be present in the 
final state, even if there were no energetic partons in the 
initial state. 

Covariant parton transport theory is the incoherent 
particle limit of the underlying nonequilibrium quantum 
theory. Its main advantages are that i) it treats both 
the low-momentum bulk sector and high-momentum rare 
particles in a self-consistent framework; ii) it is not lim¬ 
ited to local equilibrium; and iii) it describes the break-up 
of the system self-consistently. 

We consider here an inelastic extension M of the 
transport equations in Refs. 313 0 El , in which 
the on-shell parton phase space densities {fi(x,p}} evolve 
with elastic 2 —> 2 and inelastic 2 —> 2 rates as 
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\M\ 2 is the polarization averaged scattering matrix el¬ 
ement squared, the integrals are shorthands for J = 

a 

J d 3 p a /{2E a ), gi is the number of internal degrees of 
freedom for species i, while f a< i = fi(x,p a ). The source 
functions {Si(x,p)} specify the initial conditions. 

To study Au + Au collisions at RHIC at yj snn = 200 
GeV with impact parameter b = 8 fm (ss 30% central), 
we apply 0 to a system of massless gluons (g = 16), 
massless light (u,d) and strange quarks/antiquarks, and 
charm quarks/antiquarks (g = 6) with mass M c = 1.2 
GeV. The initial conditions and matrix elements were 
the same as in ID- To model the strongly-interacting 
parton system at RHIC Q, parton densities were in¬ 


creased two-fold and parton cross sections were increased 
about three-fold to dN(b=0)/drj = 2000 at midrapidity 
and Ggg^gg = 10 mb, from their perturbative estimates 
dN/dg « 1000 and a gg ^ gg « 3 mb. All elastic and in¬ 
elastic leading-order 2 —> 2 QCD processes were taken 
into account and were assumed to be enhanced by the 
same factor 10/3. The initial momentum distributions 
were taken from leading-order perturbative QCD (with 
a A'-factor of 2, GRV98LO PDFs, and Q 2 =p\, while 
Q 2 = s for charm), and the low -pr divergence in the 
light-parton jet cross sections regulated via a smooth ex¬ 
trapolation below p_i_ < 2 GeV. The transverse density 
distribution was proportional to the binary collision dis- 
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tribution for two Woods-Saxon distributions, therefore 
dN parton (b =8 frn )/dy « 500. Perfect rj = y correlation 
was assumed. 

The transport solutions were obtained via Molnar’s 
Parton Cascade algorithm m (MPC), which employs 
the parton subdivision technique 1171 to maintain Lorentz 
covariance and causality. 

Results. - There are three general ways a particle 
can acquire its final transverse momentum during the 
transport evolution: 

i) it has no interactions at all, in which case the particle 
comes from the corona of the collision region; 

ii) it loses energy via interactions and moves to lower px, 
which is the jet quenching component; or 

iii) it interacts and gains energy, i.e., gets pushed to 
higher pt by the medium. 

Figure 1 shows the relative fractions of these three 
components as a function of the final pt from the trans¬ 
port calculation for partons at midrapidity (\y/\ < 1). At 
low final pt ~ 1 — 2 GeV, partons predominantly arrive 
via the collective “push” from the medium. Contribu¬ 
tions from quenched partons are smaller, and the corona 
is negligible. On the other hand, as pr increases, the 
quenched component and the corona become more and 
more significant, while the accelerated part gradually de¬ 
creases. However, even at px ~ 7 — 8 GeV, energy loss 
gives only one half of the yield, while almost one-third of 
the partons are accelerated (initially) lower-momentum 
partons. 


come from px.i ~ 1 GeV, which is precisely where the 
initial parton px distribution (dashed line) has a strong 
peak. At the opacities considered here, partons have 
( Ncoii) ~ 14 scatterings on average. Though the like¬ 
lihood of gaining significant px , even in so many scatter¬ 
ings, is small, the initial high yield of low-momentum par- 
tons multiplied with the low probability can still compete 
with the much lower parton yields at high px . As the final 
px increases, the relative importance of low-momentum 
partons decreases, in accordance with the decrease in the 
fractional “push” yield in Fig. 1. 
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FIG. 2: Normalized initial px distributions for partons in 
various final pr intervals (and final |y/| < 1), omitting contri¬ 
butions from the corona, from the same calculation as Fig. 1. 
The initial parton pr distribution is also shown (dashed line). 
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FIG. 1: Fraction of the final high-pr yield coming from par- 
tons without interactions (corona), from energy loss (quench), 
and from accelerated partons (push) as a function of pr in 
Au + Au at y/SNN = 200 GeV with b = 8 fm, for « 6 times 
higher opacities than perturbative estimates. Computed via 

mpc BOii . 


Figure 2 gives a more detailed picture, where the nor¬ 
malized initial px distributions are shown for partons in 
various final px intervals (and final rapidity |j//| < 1), 
with contributions from the corona omitted. For all final 
px intervals, the accelerated soft partons predominantly 


Two-component models based on a thermal (or hy¬ 
drodynamic) collective part and a perturbative QCD 
jet part underestimate the collective component at high 
px because the approximately exponential px distribu¬ 
tion imposed (with typical few-hundred-MeV slope pa¬ 
rameters) cannot compete with the power-law jet spec¬ 
trum. However, such a characterization of the collective 
part is based on macroscopic approaches that are formu¬ 
lated in terms of the lowest few moments (e.g., particle 
density, energy density) of the local momentum distri¬ 
bution, and therefore cannot be trusted for the high- 
momentum tails. Our results in Figs. 1 and 2 imply 
that (initially) low-momentum partons that suffer finite 
number of scatterings populate high transverse momenta 
with a non-exponential, power-law-like distribution. In 
fact, already a single scattering generates power-law tails 
because for typical Debye-screened t-channel gluon ex¬ 
change da/dt ~ a 2 /(t — p? D ) 2 ~ a 2 /{q 2 + p? D ) 2 , where q 
is the transferred momentum. 

The above findings have significant implications for el¬ 
liptic flow. Figure 3 shows elliptic flow as a function of px 
for the three components (corona, quench, push) and for 
their combined yield. For partons that have lost energy 
(quench), elliptic flow drops rather sharply at high px, 
very similarly to calculations based on inelastic pertur- 
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bative QCD energy loss in the Eikonal limit 0 In the 
transport calculation, however, the sharp drop is largely 
compensated by the much larger V 2 of the accelerated 
partons (push). The significant “push” component also 
means that the “geometric” elliptic flow bounds argued 
in 0, which assumed extreme quenching, do not ap¬ 
ply. Though the geometric corona alone gives a positive, 
PT-independent elliptic flow (reflecting the initial spatial 
anisotropy of the edges of the collision zone), it is much 
smaller than the contribution coming from the acceler¬ 
ated partons. Therefore, the collective dynamics of the 
opaque quark-gluon plasma can provide a natural expla¬ 
nation for why elliptic flow at RHIC exceeds the “geo¬ 
metric” bounds and why it decreases relatively slowly at 
high transverse momenta [2tl| . 
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FIG. 3: Elliptic flow as a function of pr for partons without 
interactions (corona), from energy loss (quench), from acceler¬ 
ated partons (push), and from all three components combined 
(total), from the same calculation as Fig. 1. 

Conclusions. - We show from covariant parton trans¬ 
port theory that in an expanding opaque quark-gluon 
plasma there is a considerable probability that low- 
momentum partons get accelerated to high transverse 
momenta in multiple scatterings. This new mechanism of 
high-pT production, which is the opposite of jet quench¬ 
ing, has significant implications for observables at RHIC 
energies and beyond. For example, it gives a natural ex¬ 
planation to why elli ptic flow at RHIC exceeds so called 
“geometric” bounds [l9| and decreases relatively slowly 
at high transverse momenta. 

This study considered six times higher opacities than 
that of a weakly-coupled perturbative parton plasma, 
which is 2-3 times smaller than the opacities estimated 
from RHIC data (7j and likely even a smaller fraction 


of the opacities reachable at the future Large Hadron 
Collider (LHC). At larger opacities, accelerated partons 
would have even larger importance and would influence 
yields out to higher pr > 7 — 8 GeV. In addition, this 
work investigated only 2 —> 2 elastic and inelastic trans¬ 
port. A more complete study will also have to consider 
in the future the effect of radiative 2 <-> 3 processes. 
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